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I. INTRODUCTION 



The perturbative calculation of the strong nucleon-nucleon (NN) potential in the frame- 
work of effective- field theory (EFT) was first suggested by Weinberg Counting rules 
systematically arrange the magnitude of a two-nucleon irreducible diagram for calculating 
the NN potential in terms of Q/A^ where Q denotes a typical small momentum and/or pion 
mass and the chiral scale. Accordingly, one-pion exchange (OPE) and a constant two- 
nucleon contact term are the most dominant contributions. At the next-to-next-to-leading 
order (NNLO), there are two-pion exchange (TPE), contact terms with two derivatives 
and/or pion mass factors, relativistic corrections, etc. If one goes to higher orders, there 
appear multi-pion exchanges and contact terms with more than two derivatives and/or m^r, 
and heavy-meson exchange possibly comes into play. At low energies where EFT is applica- 
ble, a physical process is mostly governed by the long-range interaction and it is more or less 
insensitive to what's happening in the short-range region. Thus in most cases, heavy mesons 
are integrated out from the Lagrangian and their short-range interaction is accounted for by 
contact terms. A physical observable calculated with the A^A^ potential must be independent 
of a cutoff value that is introduced in calculating loop diagrams and its transformation to 
coordinate space. This can be achieved by contact terms together with some renormalization 
method. StrongEFT potentials thus obtained were applied to the analysis of A^A^ scattering 
phase shifts (2|, y, showing a well-behaved convergence and predictive power. The role of 
the NNLO potential was intensively explored in [sl where it was shown to be important in 
extending the predictability of the EFT to higher energies. In addition, it gives a correction 
to the leading order (LO) potential which is non-negligible even at the energy of a few MeV. 

The above approach has recently been extended to the parity-violating (PV) nucleon- 
nucleon (NN) potential [6]. Since then, it has been used for the calculation of observables 
[7, 8], especially the PV photon asymmetry in np ^7 at threshold, where it could be 
of some relevance. 

The PV A^A^ potential is obtained by replacing a parity-conserving (PC) vertex in the 
strong A^A^ potential with a PV vertex. Most of the low-energy PV calculations have relied 
on a one-meson-exchange potential with DDH estimates for the PV meson-nucleon coupling 
constants 



mo 



Theoretical estimations of have been extensively worked out with this 
del (see Refs. Q, y, Q, Q, Q for recent ones). The most elaborate results with 
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various strong-interaction models [ll|, [jj, Q, Q turn out to be basically identical. The 
asymmetry is dominated by the PV one-pion-exchange potential (OPEP) and the heavy- 
meson contribution is negligible. One can thus discuss whether the measurements of 
could provide us with an opportunity to determine the weak pion-nucleon coupling constant 
hj . Some literature has also investigated the PV two-pion-exchange potential (TPEP) 15|, 
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17l |. Its contribution to A^, calculated with the Hamada- Johnston potential, amounts to 



about —7% of the OPEP contribution 



18j. Not surprisingly, the TPEP is also part of the 



pionful EFT approach where it appears at NNLO [6||. In the present work, we concentrate 
on its contribution to the asymmetry A^ mentioned above and consider the questions that 
its estimate raises. 

Since the PV asymmetry in np d'j is sensitive to the one-pion-exchange contribution, 
we adopt the heavy-baryon chiral perturbation (HBxPT). We employ the Argonne vl8 
potential for the PC potential and the Siegert's theorem for the current operators and 
consider the PV potential relevant to "np d'y up to NNLO. This calculation will allow us 
to estimate the order and the magnitude of higher-order corrections, which will be important 
in pinning down the value of and its uncertainty. At the same time, it will provide a 
criterion for the validity of the EFT approach to the PV phenomena. 



II. FORMALISM 



In HBxPT, the order of a diagram is counted in terms of a small momentum Q (more 
precisely Q/A^) with the following rules; i) a meson propagator is counted as Q~^, ii) a 
nucleon propagator as Q~^, iii) a loop integral as Q'^, iv) a vertex as Q'^ where d is the 
number of derivatives at the given vertex. With these counting rules, diagrams in Fig. [1] 
constitute the contributions up to NNLO. One-pion exchange (a), and the two-pion exchange 
(b-d) and the contact (CT) (e) terms respectively represent LO and NNLO contributions. 

The PV potential is odd in powers of the momentum transfer in momentum space or 
the radial vector in the configuration space, which changes the orbital angular state by an 
odd number. In order to keep the wave function of a fermion system antisymmetric, the 
PV potential must satisfy the condition A{L + S + T) = even, where L, S and T denote 
orbital, spin and isospin states, respectively. The change of 5 + T in the PV potential should 
be an odd number. For a two-fermion system, the possible combinations of the spin and 
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(a) (b) (c) (d) (e) 

FIG. 1: Diagrams for PV NN potential up to NNLO: (a) for LO {0{Q-^)) OPE, and (b-e) 
for NNLO {0{Q^)) TPE + CT. Lines (dashed lines) denote nucleons (pions), vertices with a dot 
represent PC vertices, vertices with "x" represent the PV vertex proportional to h^, and a vertex 
with "®" denotes the A^A^ contact vertex function proportional to the coefficient Cg^. 

isospin operators in the PV potential must thus be {AS, AT) = (1,0) or (0,1). Among 
various possible combinations of spin and isospin operators that give A{S + T) = 1, the 
term relevant to our estimation of has the following form in momentum space 

Vi{q) = i{Ti X T2)^(o-i + 0-2) -q Vi{q), (1) 

with g = |q| and q = Pi — p'l, where pi {p[) is the momentum of a nucleon 1 in the initial 
(final) state. The OPE, TPE and CT terms are obtained as 



with 




where is the ax.al couphng constant. /. the p.ou-deeay constant and A, = 4.^ The 

constant Cq is the renormalized LEG for a A^A^ contact term Cg in |6|]. It subsumes the role 
of all the heavy degrees of freedom integrated out from the theory. We will discuss how to 
treat the renormalized LEG in our investigation after giving the formulae for the potentials 
in configuration space. 

The potential in Eq. ([1]) transformed to configuration space takes the form 
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i{ri X T2y{(Ti + 0-2) ■ [p, Vi{r)] 



(7) 



where p is the conjugate momentum of the relative coordinate r = ri — r2- For an easier 
numerical calculation, we rewrite Eqs. (0, E]) in the dispersion-relation form as 



H{q) 



°° dt' 



Ami dt' 



t'-Aml 



1 



1 1 



(8) 
(9) 



2 J 4ml Vt^ - Ami + 

In calculating the Fourier transformation of Eq. ([1]), in order to obtain a convergent, analyti- 
cal result, we introduce monopole form factors of the type (A^ — m^)/(A^ + g^) in Eq. ([2]) and 
A^/(A^ + g^) in Eqs. Hj). This particular choice for the OPEP preserves the long-range 
part, which, otherwise, could not be easily corrected for with contact interactions. The roles 
of the form factor and the cutoff are (i) to make the numerical calculation easier and more 
efficient, and (ii) to cut away the high- momentum region where the dynamics is essentially 
unknown and whose detail is irrelevant to low-energy processes. We emphasize that the 
cutoff is arbitrary for a part and that final results should not depend on its value. With the 
form factor, we rewrite the potential in configuration space as 



V/'ir) = z{n X T^na, + 0-2) ■ [p, vtir) 



where 



'-^'^ - 2V2U Anr' 



H\r)\ } , 



-Ar 



Attt ' 



with 



mlK^ r°° dt' 




A2 -t' 

g-v^r_g-Ar 



-Ar 



t' - Ami 



(10) 

(11) 
(12) 
(13) 

(14) 
(15) 



2nr J 4ml Vf ^t'-Aml ^^-i' 

Notice that in the K ^ 00 limit, L^{r) is proportional to the difference of two singular 
terms, and 5{r) with an infinite coefficient. The 3-dimensional integral over r is however 
finite and has the sign of the contact term. 



The form of depends on the regularization scheme. In Ref. Q], all the constant terms 
obtained from the dimensional regularization (d = 4 — 2e) of TPE diagrams are included in 
Cq, leading to 



C^{MX)=C,-hl^^{l-^g\) 



X 



- - 7 + ln(47r) + 2 In ( ] + 2 
e Km ' 



(16) 



where the abbreviation MX stands for 'maximal subtraction', 7 = 0.5772 and fi is the 
renormalization scale. In the minimal subtraction (MN) scheme, the renormalized LEG 
satisfies the relation 

C^iMX) = C^{MN) + C^, , (17) 

where 



c„„=-/4^(i-3ji)[ium+i 



(18) 



Either Cq{MX) or C|^(MN) has to be determined from a calculation using the underlying 
theory or from experimental data with good statistics, neither of which is available at present. 
We can only consider possible contributions to Cq, leaving for the future the theoretical or 
experimental determination of the remaining part. 

A first contribution is suggested by the one-meson exchange model. In the heavy-meson 
limit, the leading order of the p- meson propagator is a constant and the corresponding term 
can be regarded as part of Cq. It reads 



where /i^' is the PV pNN coupling constant in the potential. A non-vanishing value is ob- 
tained for this coupling in a soliton model, /i^' = —2.2 x 10"'^ [19|. Leaving aside consistency 
problems between different approaches, double counting with the TPEP considered here for 
instance, we notice that it compares to the DDH best value of /i^ (= 4.6 x 10^^). The 
corresponding value of C^{p) is — L20 x lO^^MeV^'^ in units of this value of h^.. 

A second contribution is provided by the MN scheme where Cmn is treated independently 
of Cff (MN), which consequently gives p dependence in the result. If we choose /i = nip, we 
have Cmn = 4.174 x 10~^MeV~^ in units of /i^. This value is one order larger than C|^(p), 
and therefore the effect of Cmn and the dependence on p can be non-negligible. 

A third contribution stems from a deeper examination of the TPEP. This one has been 
obtained by removing from the square-box diagram the iterated OPEP, ignoring the cutoff 



introduced in actual calculations. The correction, which involves the whole strong interaction 
and should contribute to make the results cutoff independent, is not easy to calculate. This 
last property can however be restored with a minimal contribution which cancels the cutoff- 
dependent term in Eq. flTTl) . This one corresponds to 

Its value compares to the Cmn one for the smallest values of the cutoff A considered here. 

We will explore the above contributions in detail in the analysis of the photon asymmetry 
in np — > d'j, A^. This quantity is defined from the differential cross section of the process as 

^ oc 1 + y4^cos6', (21) 
ail 

where 6 is the angle between the neutron polarization and the out-going photon momentum. 
Non-zero A^ values arise from the interference of opposite-parity transition amplitudes, e.g. 
Ml and El. At the thermal energy where the process occurs, lowest order EM operators 
may suffice, therefore we consider the El operator, 

Jei = -z^K-r|)r, (22) 

where is the out-going-photon energy. At the leading order of h^, A^ is proportional to 
hi, and we can write A^ as 

A^ = a^hl , (23) 

with 

_ _ ReCM^ 

where Si and Aii are matrix elements of the El and Ml transitions, respectively. Analytic 
forms of these amplitudes can be found in |14]. 

III. RESULTS AND DISCUSSION 

Here we present and discuss numerical results corresponding to the various contributions 
outlined in the previous section. The OPE and TPE potentials in configuration space are 
considered as well as individual and total contributions to the A^ asymmetry. 
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FIG. 2: Potentials Vi^{r) (upper row) and V2T^{r) (lower row) divided by h\ in the short- 
intermediate (left column) and long (right column) ranges. In the figures for TPEP, the result 
for A = 5000 MeV is shown as an example of the limiting case A ^ oo. 

The OPEP term v^^{r) is shown in the top rows of Fig. [2]in the short (left) and long (right) 
ranges separately. The curve denoted by 'Yukawa' corresponds to the infinite-cutoff value. 
Increasing the cutoff A, the potential converges to the Yukawa one. While the potential is 
almost independent of the cutoff value in the range r > 1 fm, some sizable deviation from 
the Yukawa form occurs at r < 1 fm when the cutoff is equal to 500 MeV. 

The TPEP term f^^(r) is in the bottom row of Fig. [2J The long-range behavior is, 
as OPEP, independent of the cutoff value, but its magnitude is smaller by an order of 
magnitude. Therefore one can conclude that the long-range dynamics is essentially cutoff 
independent and governed by OPEP. On the contrary, TPEP changes sign in the interme- 
diate or short range (bottom-left panel), and the position at which the sign change occurs 
depends on the cutoff value. What matters in this behavior is that this change occurs in the 
intermediate region where TPEP is supposed to play an important role. The effect of this 
strong cutoff dependence is investigated below for the physical observable of interest here. 
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A (MeV) 500 1000 1500 2000 



„OPE 




nQQ9 


n 1 1 DA 


n 1 1 1 7 


n 1 1 1 Q 

— u. Illy 






-0.0022 


0.0073 


0.0117 


0.0133 




(P) 


-0.0008 


-0.0004 


-0.0002 


-0.0001 


„LEC 


(mn) 


0.0264 


0.0134 


0.0065 


0.0036 




(7r,A) 


-0.0128 


-0.0016 


-0.0004 


-0.0001 



TABLE I: Contribution to the asymmetry from each term as a function of the cutoff A. The 
renormahzation scale is chosen as = nip in the minimal subtraction LEG contribution (oc Cmn )• 

Table [T] shows the contribution of each term to the coefficient entering the expression 
of the asymmetry A^, Eq. (123|) . The OPEP contribution with A = 500 MeV is slightly 
smaller in magnitude than those with the other cutoff values, but the result as a whole 
is independent of the cutoff value. The small dependence of a°^^ on the cutoff A can be 
attributed to the sensitivity of v^^^r) in the range r < 1 fm to A. On the other hand, it 
seems that the region below r ~ 0.5 fm gives an almost negligible contribution since Vi^{r) 
with A = 1000 MeV differs sizably from its value with larger cutoff values below r ~ 0.5 fm, 
but the difference in a^^^ is negligible. Roughly speaking, the significant contribution to 
comes from the region r ^ 0.5 fm. 

The TPEP contribution to a^, on the other hand, varies significantly in sign and magni- 
tude. The change of a^^^ can be attributed to the strong dependence of f^^(r) on the cutoff 
value in the intermediate and short ranges. The sign change, which can be ascribed to a 
contact term implied by the method used to calculate the potential, is physically question- 
able. From the behavior of v^^^r) and the sign of a^^^, one can deduce that a decreasing 
term gives a negative contribution to a^. With a smaller cutoff value, the decreasing part 
of TPEP spreads to a larger r. For instance, TPEP becomes increasing at around r = 0.4 
fm for A = 2000 MeV, but this occurs around 1.5 fm for A = 500 MeV. In the discussion of 
a^^^, we argued that most of the contribution arises from the region r ^ 0.5 fm. With a 
smaller cutoff value, the negative contribution to a^^^ becomes more substantial and, as a 
result, aj^^^ with A = 500 MeV is negative. The convergence of the TPE contribution with 
an increasing cutoff value is slower than the OPE one, which can also be attributed to the 
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A (MeV) 500 1000 1500 2000 



a*°*(MX) -0.1149 -0.1051 -0.1006 -0.0988 

H = 2m^ -0.0987 -0.0969 -0.0966 -0.0966 

a^°*(MN) /u = mp -0.0885 -0.0917 -0.0941 -0.0952 
H = l GeV -0.0859 -0.0904 -0.0935 -0.0949 



TABLE II: Total sum of OPE, TPE and LEG contributions in the MX and MN schemes. Each 
row of aIj°*(MN) corresponds to fj, = 2m,r> fnp and 1000 MeV from the top, respectively. 

strong dependence of TPEP on the cutoff value. However, a!^^^ will not increase indefinitely 
and we guess that a^^^ with A = 2000 MeV will not be much different from that in the 
limiting case A oo. For example, aj^^^ = 0.0146 with A = 5000 MeV, which suggests 
that the rate of increase with respect to A has decreased compared to that in the smaller A 
region. Thus, without significant error, we can conclude that the TPE contribution to is 
about 10% of the OPE. 

All the LEG contributions (oc Cq{p), Cmn and C(f (vr. A)) have a magnitude that decreases 
with increasing cutoff value. With a smaller cutoff value, Vq^j-^v) is less short-range peaked 
and more spread out to larger r. Taking into account the r ^ 0.5 fm criterion, it is a natural 
result to have a a^^^ smaller in magnitude with a larger cutoff. The Cq{p) contribution 
is small in comparison of the other ones, but this could be a result of the particular value 



of h^p used in its calculation [19|. The contribution of the Cmn term, absent in the MX 
scheme, is comparable to the TPEP one. Consequently, the choice of the regularization 
scheme can give a non-trivial effect for the total sum of OPE, TPE and LEG contributions, 
al°^. The Cq{tt,A) contribution is also comparable to the TPEP one, but only at small A. 
As expected from its derivation, it removes the A dependence of the OPE contribution. 

In Table HIl we show the sum of OPE, TPE and LEG contributions to a^, for the MX 
and MN regularization schemes. Comparison of the results in different columns gives insight 
on the cutoff dependence, while comparion of rows in the MN case shows how results vary 
with the renormalization scale, fi. This scale is conventionally set to the p-meson mass. 



the chiral symmetry breaking scale A^ or the order of the cutoff value A Exploring 
the corresponding dependence with values fi = 2m.,r, rrip and 1000 MeV, a!j°* is found to 
increase with p (algebraically). This is simply related to the contribution of the Cmn term 
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which is a logarithmically increasing function of that variable. The dependence on the cutoff 
value is found to be opposite in the MX scheme and the MN one with fj, — nip. Since the 
contributions of the OPE plus Cq{tt, A) and C^{p) terms remain almost constant, the A 
dependence in the MX scheme is mainly due to the TPE one. The opposite effect in the 
MN scheme results from the further LEG contribution involving Cmn, which roughly has the 
same magnitude but is decreasing instead of increasing with A. The overall A dependence 
is sensitive to the ^ value and can almost vanish (// = 2m7r) as well as be more pronounced 
1000 MeV). 

Because of the considerable dependence on the cutoff A (MX scheme) or on both the 
cutoff and the renormalization scale (MN scheme), the uncertainty of a!^°*(MX) or a!^°*(MN) 
is rather large. Moreover the average value differs from one scheme to the other. In principle, 
LEG contributions should make low-energy observables independent of the input parameters. 
In this respect, our results are not quite satisfactory. In the MX scheme, further contributions 
are definitively required, while in the MN scheme, instead, it appears that an acceptable 
result could be obtained, provided that a low value of is used. Interestingly, this occurs 
for a value of Cmn that roughly compensates the TPEP contribution at short distances, 
which is partly unphysical. A cancellation on this basis implies /i = A in the log term 
contributing to Cmn (and perhaps some modification of the associated constant) . Whatever 
the approach, results show that OPE is a dominant contribution to a^, and that the NNLO 
contributions are in a reasonable range of about 10 ~ 20 % of the OPE contributions for the 
most stable cases. We would say it reasonable in the sense that the range is obtained with 
standard values of input parameters and that, at the same time, the amount of the NNLO 
contribution has the typical size of higher-order corrections in the EFT calculation. 

IV. SUMMARY 

We have considered the PV NN potential, up to the NNLO, from heavy-baryon chiral- 
perturbation theory, and applied it to the calculation of the PV asymmetry in np — > 
at threshold. The OPEP appears at the LO, there is no term at the NLO, and TPEP 

and GT terms are picked up at the NNLO. Heavy degrees of freedom are integrated out by 
introducing a monopole form factor in the Fourier transformation of the momentum-space 
potentials and the corresponding contribution is ascribed to contact terms. A renormalized 
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LEG is obtained from the dimensional regularization and its form depends on the regula- 
rization scheme. We employ the minimal and maximal subtraction schemes to determine 
this form. Lacking a way to fix the renormalized LEG, we illustrated its contribution by 
two terms. The first one is obtained from the p-resonance saturation of the vector-meson 
exchange PV potential, and the other one restores the A independence of the OPE con- 
tribution. In the MN scheme, we have an additional constant term which depends on the 
renormalization scale. Gutoff dependence is explored with values A = 500 ~ 2000 MeV, and 
regularization-scheme dependence with renormalization scales /i = 2m^ ~ 1000 MeV. 

The OPE contribution to the asymmetry is satisfactorily cutoff independent. The TPE 
potential depends strongly on the cutoff value and is somewhat uncertain, in particular 
in the intermediate range. As a result, the TPE contribution to the asymmetry varies 
widely and even shows a change of sign. However, the TPE contribution is bounded by the 
limiting cutoff value A oo and does not exceed about 10% of the OPE one. The LEG 
contributions are found to be significantly dependent on both the cutoff and renormalization 
scales introduced in their calculations. 

The total sum of OPE, TPE and LEG contributions varies non-negligibly, depending on 
the choice of the cutoff and renormalization-scale values. With the parameters considered, 
we have alj°* = —(0.08 ~ 0.11). If the PV "np d'j experiment measures at the order 
10~^, the present result of al^°* can fix the first digit of the weak tcNN coupling constant /i^ 
unambiguously. However, some caution is needed. In principle, the role of LEG is to account 
for degrees of freedom whose microscopic description is practically irrelevant for low-energy 
processes like the one considered here and, at the same time, to make results independent of 
parameters or regularization schemes employed in their calculation. The present treatment 
of the LEG doesn't fully satisfy the latter requirement, though some steps in this direction 
have been considered. Thus, the estimate of al^°* made here provides a range in which its 
true value may be, but it is still premature to claim the result for sure. For pinning down 
the true value of alj°*, the LEG terms should be treated in a more rigorous and reliable way, 
e.g. fitting them to experimental data if possible, obtaining them with lattice calculations, 
or making a better use of renormalization group methods. A comparison with earlier TPEP 



calculations 
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17| could also be quite useful. 
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Abstract 

The parity-violating nucleon-nucleon {NN) potential is considered up to next-to-next-to leading 
order in heavy-baryon chiral perturbation theory. We include the one-pion exchange at the leading 
order and the two-pion exchange and two-nucleon contact terms at the next-to-next-to-leading 
order. The effects of intermediate (two-pion exchange) and short-range (two-nucleon contact) terms 
are probed by calculating the photon asymmetry in ftp — )■ dj employing Siegert's theorem and 
an accurate phenomenological potential for the parity- conserving NN interaction. We explore in 
detail the uncertainties due to the parameters that control the contribution of the short-range 
interaction. We obtain about 20% uncertainty in the value of Aj up to the next-to-next-to leading 
order. We discuss its implication for the determination of the weak pion-nucleon coupling constant 
and how the uncertainty can be reduced. 
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I. INTRODUCTION 



The perturbative calculation of the strong nucleon-nucleon (NN) potential in the frame- 
work of effective-field theory (EFT) was first suggested by Weinberg Counting rules 
systematically arrange the magnitude of a two-nucleon irreducible diagram for calculating 
the A^A^ potential in terms of Q/A^ where Q denotes a typical small momentum and/or pion 
mass and the chiral scale. Accordingly, one-pion exchange (OPE) and a constant two- 
nucleon contact term are the most dominant contributions. At the next-to-next-to-leading 
order (NNLO), there are two-pion exchange (TPE), contact terms with two derivatives 
and/or pion mass factor, relativistic corrections, etc. If one goes to higher orders, there 
appear multi-pion exchanges and contact terms with more than two derivatives and/or m^, 
and heavy-meson exchange possibly comes into play. At low energies where EFT is applica- 
ble, a physical process is mostly governed by the long-range interaction and it is more or less 
insensitive to what's happening in the short-range region. Thus in most cases, heavy mesons 
are integrated out from the Lagrangian and their short-range interaction is accounted for by 
contact terms. A physical observable calculated with the A^A^ potential must be independent 
of a cutoff value that is introduced in calculating loop diagrams and its transformation to 
coordinate space. This can be achieved by contact terms together with some renormalization 



method. Strong E 
phase shifts 



^'T potentials thus obtained were applied to the analysis of A^A^ scattering 
J], showing a well-behaved convergence and predictive power. The role 



of the NNLO potential was intensively explored in j5| where it was shown to be important 
in extending the predictability of the EFT to higher energies. In addition, its correction to 
the leading order (LO) potential is non-negligible even at the energy of a few MeV. 

The above approach has recently been extended to the parity-violating (PV) nucleon- 
nucleon (NN) potential [g]. Since then, it has been used for the calculation of observables 



a a, 



especially the PV photon asymmetry Aj in np dj at threshold, where it could be 
of some relevance. 

The PV A^A^ potential is obtained by replacing a parity-conserving (PC) vertex in the 
strong A^A^ potential with a PV vertex. Most of the low-energy PV calculations have been 
relying on a one-meson-exchange potential with DDH estimates for the PV meson-nucleon 
coupling constants j^. Theoretical estirnations of have been extensively worked out with 



n 



this model (see Refs. [lOl, mi, Q, y , Q for recent ones). The most elaborate results with 
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various strong-interaction models [ill, Q, Q, [l^ turn out to be basically identical. The 
asymmetry is dominated by the PV one-pion-exchange potential (OPEP) and the heavy- 
meson contribution is negligible. One can thus discuss whether the measurements of Aj 
could provide us with an opportunity to determine the weak pion-nucleon coupling constant 
hi. Some literature investigated the PV two-pion-exchange potential (TPEP) in the past 



16 



13]. Its contribution to A^, calculated with the Hamada- Johnston potential, amounts to 



about —7% of the OPEP contribution [18|. Not surprisingly, the TPEP is also part of the 
pionful EFT approach where it appears at NNLO js]. In the present work, we concentrate 
on its contribution to the asymmetry A^ mentioned above and consider the questions that 
its estimate raises. 

Since the PV asymmetry in np dj is sensitive to the one-pion-exchange contribution, 
we adopt the heavy-baryon chiral perturbation (HBxPT). We employ the Argonne vl8 
potential for the PC potential and the Siegert's theorem for the current operators and 
consider the PV potential relevant to np d'j up to NNLO. This calculation will allow us 
to estimate the order and the magnitude of higher-order corrections, which will be important 
in pinning down the value of and its uncertainty. At the same time, it will provide a 
criterion for the validity of the EFT approach to the PV phenomena. 



II. FORMALISM 



In the HBxPT, the order of a diagram is counted in terms of a small momentum Q (more 
precisely Q/A^) with the following rules; i) a meson propagator is counted as Q^^, ii) a 
nucleon propagator as Q^^, iii) a loop integral as Q"^, iv) a vertex as Q'^ where d is the 
number of derivatives at the given vertex. With these counting rules, diagrams in Fig. [T] 
constitute the contributions up to NNLO. One-pion exchange (a), and the two-pion exchange 
(b-d) and the contact (CT) (e) terms respectively represent LO and NNLO contributions. 

The PV potential is odd in powers of the momentum transfer in momentum space or 
radial vector in the configuration space, which changes the orbital angular state by an odd 
number. In order to keep the wave function of a fermion system antisymmetric, the PV 
potential must satisfy the condition A{L + S + T) = even, where L, S and T denote orbital, 
spin and isospin states, respectively. The change of S + T in the PV potential should be an 
odd number. For a two-fermion system, the possible combinations of the spin and isospin 
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(a) (b) (c) (d) (e) 

FIG. 1: Diagrams for PV NN potential up to NNLO: (a) for LO (C'(Q"^)) OPE, and (b-e) 
for NNLO {0{Q^)) TPE + CT. Lines (dashed lines) denote nucleons (pions), vertices with a dot 
represent PC vertices, vertices with "x" represent the PV vertex proportional to h^, and a vertex 
with "(E)" denotes the NN contact vertex function proportional to the coefficient Cff . 

operators in the PV potential must thus be {AS, AT) = (1,0) or (0,1). Among various 
possible combinations of spin and isospin operators that give A{S + T) = 1, the term 
relevant to our estimation of has the following form in momentum space 

Vi{q) = i{Ti X T2y{(Ti + 0-2) -q Vi{q), (1) 

where q = \q\ and q = Pi — P2- The OPE, TPE and CT terms are obtained as 



with 




where is the axial coupling constant, /. the pion-deoay constant and A, = 4./,. The 

constant is the renormalized LEG for a NN contact term Ce in p. It subsumes the role 
of all the heavy degrees of freedom integrated out from the theory. We will discuss how to 
treat the renormalized LEG in our investigation after giving the formulae of the potentials 
in configuration space. 

The potential in Eq. ([T]) transformed to the configuration space takes the form 

= z(ti X T2Y{(Ti + (T2) • [p, Vi{r)] , (7) 

4 



where p is the conjugate momentum of the relative coordinate r = ri — r2. For an easier 
numerical calculation, we rewrite Eqs. (ISl [6|) into the dispersion-relation form as 

°° dt' 



~ 4ml f°° dt' 1 1 

^^^^ = Li Tv^t'-Amlt' + g' ■ 

In calculating the Fourier transformation of Eq. and in order to make it analytically, we 
introduce monopole form factors of the type (A^ — m^) / (A^ + g^) in Eq. ^ and A^/ (A^ + q"^) 
in Eqs. ^ [4|. The particular choice for the OPEP preserves the long-range part, which, 
otherwise, could not be easily corrected for with contact interactions. The roles of the form 
factor and the cutoff are (i) to make the numerical calculation easier and more efficient, and 
(ii) to cut away the high-momentum region where the dynamics is essentially unknown and 
whose detail is irrelevant to low-energy processes. It is reminded that the cutoff is arbitrary 
for a part and that final results should not depend on its value. With the form factor, we 
rewrite the potential in configuration space as 



where 



with 



V/'ir) = z{n X r,y{a^ + a^) ■ [p, vt{r)] , (10) 



vtir) = v^TT^ {g^L^r) - g\ [ZL\r) - H\r)\ } , (12) 

X 

— Kt 

v^Ar) = -C^^'"-^ . (13) 



2a2 „oo fit' 1 p-v^^'—p-Ar 

^ ' 27rr J,mi v^./^riw ^'-t' ' ^ ' 



Notice that in the A ^ oo limit, L^{r) is the difference of terms as singular as r^^ and 5{r) 
with an infinite coefficient. The 3-dimensional integral over r is however finite and has the 
sign of the contact term. 
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The form of Cjf depends on the regularization scheme. In Ref. j^, all the constant terms 
obtained from the dimensional regularization (d = 4 — 2e) of TPE diagrams are included in 
C^, leading to 

C^iMX) = Ce- - 3g\) - 7 + ln(47r) + 2 In (-^) + 2] , (16) 

where the abbreviation MX stands for 'maximal subtraction', 7 = 0.5772 and ji is the 
renormalization scale. In the minimal subtraction (MN in short), the renormalized LEG 
satisfies the relation 

C^iMX) = C^iMN) + C™„ , (17) 

where 



(18) 



Either C^(MX) or Cjf (MN) has to be determined from the calculation with the underlying 
theory or experimental data with good statistics but none of them is available at present. 
We can only consider possible contributions to C^, leaving for the future the theoretical or 
experimental determination of the remaining part. 

A first contribution is suggested by the one-meson exchange model. In the heavy-meson 
limit, the leading order of the p-meson propagator is a constant and the corresponding term 
can be regarded as part of C^. It reads 

^«<'" = f^- (1^) 

where hjl is a PV pNN coupling constant in the potential. A non- vanishing value is obtained 
for this coupling in a soliton model, /i^' = —2.2 x 10~^ [l^. Leaving aside consistency 
problems between different approaches, double counting with TPEP considered here for 
instance, we notice that it compares to the DDH best value of (= 4.6 x 10^^). The 
corresponding value of Cq{p) is —1.20 x 10^^ MeV^^ in unit of this value of h]^. 

A second contribution is provided by the MN scheme where Cmn is treated independently 
of C^(MN), which consequently gives fi dependence in the result. If we choose // = rUp, we 
have = 4.174 x 10"^ MeV"^ in unit of h^. This value is one order larger than C^{p), 
and therefore the eff'ect of Cmn and the dependence on ji can be non-negligible. 

A third contribution stems from a deeper examination of the TPEP. This one has been 
obtained by removing from the square-box diagram the iterated OPEP, ignoring the cutoff 

6 



introduced in actual calculations. The correction, which involves the whole strong interaction 
and should contribute to make the results cutoff independent, is not easy to calculate. This 
last property can however be restored with a minimal contribution which cancels the cutoff- 
dependent term in Eq. (ITTi ). This one corresponds to 

^'<-^' = -5^- 

Its value compares to the one for the smallest values of the cutoff A considered here. 

We will explore the above contributions in detail in the analysis of the photon asymmetry 
in np — > d'j, A^. This one is defined from the differential cross section of the process as 

— oc 1 + cos 6*, (21) 

where 9 is the angle between the neutron polarization and the out-going photon momentum. 
Non-zero A^ values arise from the interference of opposite-parity transition amplitudes, e.g. 
Ml and El. At the thermal energy where the process occurs, lowest order EM operators 
may suffice, therefore we consider the El operator, 

Jei = (r^ - r|) r , (22) 

where cOry is the out-going-photon energy. At the leading order of h^, A^ is proportional to 
hi, and we can write Aj as 

A^ = a^hl , (23) 

with 

where £i and A4i are matrix elements of the El and Ml transitions, respectively. Analytic 
forms of these amplitudes can be found in fl^ . 

III. RESULTS AND DISCUSSION 

We here present and discuss numerical results corresponding to various contributions 
outlined in the previous section. The OPE and TPE potentials in configuration space are 
considered as well as individual and total contributions to the A^ asymmetry. 
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FIG. 2: Potentials Ui^(r) (upper row) and V2Tj{r) (lower row) divided by h\ in the short- 
intermediate (left column) and long (right column) ranges. In the figures for TPEP, the result 
of A = 5000 MeV is shown as an example of the limiting case A — >■ oo. 

The OPEP term v^^(r) is shown in the top rows of Fig. [2] in the short (left) and long 
(right) ranges separately. The curve denoted by 'Yukawa' corresponds to the infinite-cutoff 
value. Increasing the cutoff A, the potential converges to the Yukawa one. It is almost 
independent of the cutoff value in the range r > 1 fm and deviates from the others at r < 1 
fm with a cutoff of 500 MeV, and at roughly r = 0.5 fm with a cutoff of 1000 MeV. 

The TPEP term V2T^{r) is in the bottom rows of Fig. [2l The long-range behavior is, 
as OPEP, independent of the cutoff value, but its magnitude is smaller by an order of 
magnitude. Therefore one can conclude that the long-range dynamics is essentially cutoff 
independent and governed by OPEP. On the contrary, TPEP changes sign in the interme- 
diate or short range (bottom-left panel), and the position at which the sign change occurs 
depends on the cutoff value. What matters in this behavior is that this change occurs in the 
intermediate region where TPEP is supposed to play an important role. The effect of this 
strong cutoff dependence is investigated below for the physical observable of interest here. 
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A (MeV) 500 1000 1500 2000 
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n 0009 
— u.uyyz 


n 1 1 0A 

— U. i IU'4: 


1117 
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n 1 1 1 

— u. 1 i ly 




-0.0022 


0.0073 


0.0117 


0.0133 




-0.0008 


-0.0004 


-0.0002 


-0.0001 




0.0264 


0.0134 


0.0065 


0.0036 




-0.0128 


-0.0016 


-0.0004 


-0.0001 



TABLE I: Contribution to the asymmetry from each term as a function of the cutoff A. The 
renormalization scale is chosen as = mp in the minimal subtraction LEG contribution (oc Cmn )■ 

Table [T] shows the contribution of each term to the coefficient entering the expression 
of the asymmetry A^, Eq. The OPEP contribution with A = 500 MeV is slightly 

smaller in magnitude than those with the other cutoff values, but the result as a whole is 
independent of the cutoff value. The small difference of a°^^ with A = 500 MeV can be 
attributed to the deviation of Vi^((r) in the range r < 1 f m with A = 500 MeV from those 
with the other cutoff values. On the other hand, it seems that the region below r ~ 0.5 fm 
gives an almost negligible contribution since v^^((r) with A = 1000 MeV differs sizably from 
those with larger cutoff values below r ~ 0.5 fm, but the difference for a^^^ is negligible. 
Roughly speaking, the significant contribution to comes from the region r ^ 0.5 fm. 

The TPEP contribution to a^, on the other hand, varies significantly in sign and magni- 
tude. The change of d^^^ can be attributed to the strong dependence of t'2V(r) on the cutoff 
value in the intermediate and short ranges. The sign change, which can be ascribed to a 
contact term implied by the method used to calculate the potential, is physically question- 
able. ^From the behavior of v^^{r) and the sign of a°^^, one can deduce that a decreasing 
term gives a negative contribution to a^. With a smaller cutoff value, the decreasing part 
of TPEP spreads to a larger r. For instance, TPEP becomes increasing at around r = 0.4 
fm for A = 2000 MeV, but this occurs around 1.5 fm for A = 500 MeV. In the discussion of 
a°^^, we argued that most of the contribution arises from the region r ^ 0.5 fm. With a 
smaller cutoff value, the negative contribution to d^^^ becomes more substantial and, as a 
result, a^^^ with A = 500 MeV is negative. The convergence of the TPE contribution with 
an increasing cutoff value is slower than the OPE one, which can also be attributed to the 
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A (MeV) 500 1000 1500 2000 



a\"\MX) -0.1149 -0.1051 -0.1006 -0.0988 
-0.0987 -0.0969 -0.0966 -0.0966 

a^°*(MN) -0.0885 -0.0917 -0.0941 -0.0952 
-0.0859 -0.0904 -0.0935 -0.0949 



TABLE II: Total sum of OPE, TPE and LEG contributions in the MX and MN schemes. Each 
row of al^°*(MN) corresponds to ^ = 2m7r, nip and 1000 MeV from the top, respectively. 

strong dependence of TPEP on the cutoff value. However, d^^^ will not increase indefinitely 
and we guess that a!^^^ with A = 2000 MeV will not be much different from that in the 
limiting case A ^ oo. For example, aj^^^ = 0.0146 with A = 5000 MeV, whose increase rate 
with respect to A is very lowered down compared to that in the smaller A region. Thus, 
without significant error, we can conclude that the TPE contribution to is about 10% of 
the OPE. 

All the LEG contributions (oc Cff (p), Cmn and Cff (tt. A)) have a magnitude that decreases 
with increasing cutoff value. With a smaller cutoff value, VQrj,{r) is less short-range peaked 
and spreads broader to larger r. Reminding the r ^ 0.5 fm criterion, it is a natural result 
to have a a^^'^ smaller in magnitude with a larger cutoff. It is noticed that the C^{p) 
contribution is small in comparison of the other ones but this could be a result of the 
particular value of used in its calculation [l9]. The contribution of the Cmn term, absent 
in the MX scheme, compares to the TPEP one. Consequently, the choice of the regularization 
scheme can give a non-trivial effect for the total sum of OPE, TPE and LEG contributions, 
^tot fjij^g C^(7r,A) contribution also compares to the TPEP one, but only at small A. As 
expected from its derivation, it removes the A dependence of the OPE contribution. 

In Table HTl we show the sum of OPE, TPE and LEG contributions to a^, for the MX and 
MN regularization schemes. Columns give insight on the cutoff dependence while rows in the 
MN case show how results vary with the renormalization scale, //. This one is conventionally 
set to the p-meson mass, the chiral symmetry breaking scale A^ or also the order of the 
cutoff value A Q]. Exploring the corresponding dependence with values /i = 2mT^, rUp 
and 1000 MeV, al^°* is found to increase with jj, (algebraically). This is simply related 
to the contribution of the Cmn term which is a logarithmically increasing function of that 
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variable. Concerning the dependence on the cutoff value, it is found to be opposite in the 
MX scheme and the MN one with /j, = nip. As the OPE together with the C^{'k, A) and 
Cq{p) contributions remains almost constant, the A dependence in the MX scheme is mainly 
due to the TPE one. The opposite effect in the MN scheme results from the further LEG 
contribution involving Cmn, which roughly has the same magnitude but is decreasing instead 
of increasing with A. It is noticed that the overall A dependence is sensitive to the fi value 
and can almost vanish (// = 2m^) as well as be more pronounced {fi = 1000 MeV). 

Because of the considerable dependence on the cutoff A (MX scheme) or on both the 
cutoff and the renormalization scale (MN scheme), the uncertainty of a^°*(MX) or a^°*(MN) 
is rather large. Moreover the average value differs from one scheme to the other. In principle, 
LEG contributions should make low-energy observables independent of the input parameters. 
In this respect, our results are not quite satisfactory. In the MX scheme, further contributions 
are definitively required, while in the MN scheme, instead, it sounds that an acceptable result 
could be obtained, provided that a low value of fi is used. Interestingly, this occurs for a 
value of Cmn that roughly compensates the TPEP contribution at short distances, which is 
partly unphysical. A cancellation on this basis implies /i = A in the log term contributing to 
Cjnn (and perhaps some modification of the associated constant). Whatever the approach, 
results show that OPE is a dominant contribution to a^, and that the NNLO contributions 
are in a reasonable range of about 10 ~ 20 % of the OPE one for the most stable cases. 
We would say it reasonable in the sense that the range is obtained with standard values of 
input parameters and that, at the same time, the amount of the NNLO contribution has 
the typical size of higher-order corrections in the EFT calculation. 

IV. SUMMARY 

We have considered the PV NN potential up to the NNLO from the heavy-baryon chiral- 
perturbation theory, and applied it to the calculation of the PV asymmetry in ftp at 
threshold. The OPEP appears at the LO, there is no term at the NLO, and TPEP and 
CT terms are picked up at the NNLO. Heavy degrees of freedom are integrated out by 
introducing a monopole form factor in the Fourier transformation of the momentum-space 
potentials and the corresponding contribution is ascribed to contact terms. A renormalized 
LEG is obtained from the dimensional regularization and its form depends on the regulari- 
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zation scheme. We employ the minimal and maximal subtraction schemes to determine this 
form. Lacking a way to fix the renormalized LEG, we illustrated its contribution, for a part, 
with a constant obtained from the p-resonance saturation of the vector-meson exchange PV 
potential. Another part restores the A independence of the OPE contribution. In the MN 
scheme, we have an additional constant term which depends on the renormalization scale. 
Cutoff dependence is explored with values A = 500 ~ 2000 MeV, and regularization-scheme 
dependence with renormalization scales /i = 2771,^ ~ 1000 MeV. 

The OPE contribution to the asymmetry is satisfactorily cutoff independent. The TPE 
potential depends strongly on the cutoff value and is somewhat uncertain, in particular 
in the intermediate range. As a result, the TPE contribution to the asymmetry varies 
widely and even shows a change of sign. However, the TPE contribution is bounded by the 
limiting cutoff value A ^ oo and does not exceed about 10% of the OPE one. The LEG 
contributions are found to be significantly dependent on both the cutoff and renormalization 
scale introduced in their calculations. 

The total sum of OPE, TPE and LEG contributions varies non-negligibly, depending on 
the choice of the cutoff and renormalization-scale values. With the considered parameters, 
we have a^°* = —(0.08 ~ 0.11). If the PV np experiment measures at the order 

10^^, the present result of a^°* can fix the first digit of the weak ttNN coupling constant 
unambiguously. However, some caution is needed. In principle, the role of LEG is to account 
for degrees of freedom whose microscopic description is practically irrelevant for low-energy 
processes like the one considered here and, at the same time, to make results independent of 
parameters or regularization schemes employed in their calculation. The present treatment 
of the LEG doesn't fully satisfy the latter requirement, though some steps in this direction 
have been considered. Thus, the estimate of a^°* made here provides a range in which its 
true value may be, but it is still premature to claim the result for sure. For pinning down 
the true value of al^°*, the LEG terms should be treated in a more rigorous and reliable way, 
e.g. fitting them to experimental data if possible, obtaining them with lattice calculations, 
or making a better use of renormalization group methods. A comparison with earlier TPEP 



calculations 
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13] could also be quite useful. 
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